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I. INTRODUCTION 



It is a commonly accepted scenario for central ultra-relativistic heavy-ion collisions that immediately after the 
collision, sandwiched between the receding baryonic pancakes, a rapidly expanding spatial region is left, characterized 
by high energy density and low baryon density 0. The expansion may be mainly in beam direction, such that 
the longitudinal velocities of volume elements are proportional to their distance from the collision point. Assuming 
invariance with respect to transformations between comoving frames, the rapidity density of particles emitted from 
this expanding rod will show a plateau reflecting that symmetry. This plateau is experimentally well established, and 
recent results for baryonic yields indicate that for high centrality and increasing energy the antibaryon/baryon ratio 
in fact tends towards unity [3- 0- C3| ■ 

The production of baryon- antibaryon pairs from regions of highly excited 'vacuum' presents a challenge for many 
different types of theoretical models. It arises very naturally in models where baryon density is identified with 
topological winding density in effective meson field theories 0> 0, S ■ F° r random field configurations on random 
triangular lattices the probability for nontrivial winding provides an upper (Kibble) limit for baryon-antibaryon 
multiplicities [3. Ilfjj. The dynamical evolution in time, which follows the preparation of a statistical ensemble at 
some initial time, tends to smoothen the random initial configurations, resulting in frequent 'annihilation' of positive 
and negative winding density, such that the final baryon-antibaryon multiplicities will strongly depend on the cooling 
mechanism and freeze-out ti mes, after which the system is sufficiently decoupled, such that baryons and antibaryons 
no longer interact [TlLIT3.ll3.ll4j. 

Apart from genuine energy dissipation through meson emission, the rapid expansion provides a kinetic cooling 
mechanism which naturally decreases all gradients in a field configuration and, consequently, the winding densities. 
This expansion is conveniently accounted for by transforming the space-time frame into the ra pidi ty - proper-time 
frame and formulating the field dynamics in these locally comoving frames [TEL 

EmEmMiiTi. We repeat 

the essential points of this formulation in section II, with specific attention to the soliton solutions. We also include 
the effects of genuine damping and discuss the adiabatic approximation which applies for overdamped systems and 
provides significant simplifications in the numerical simulation of such evolutions. 

In section III of the present paper we apply this formalism to the (1+1) dimensional <i> 4 model, which is characterized 
by only two discrete degenerate vacua. This provides a very transparent example for the mechanisms which govern the 
'phase transition' from the hot initial configurations through the freeze-out into the final cold system which has settled 
into the true vacua and consists of non-interacting solitons and vacuum fluctuations. We obtain simple analytical 
expressions for freeze-out time and kink multiplicities, both for the system moving freely in the expanding Bjorken 
rod, and for the damped system which additionally looses energy through a dissipative term. 

II. FIELD DYNAMICS IN THE RAPIDITY - PROPER-TIME FRAME 
A. Transformation to Bjorken frame 

As standard effective action S for a scalar field $ in space-time (z, t) we consider 

S = J Q(d t $) 2 - \{d z ^f - V(*)j dzdt (1) 

where V(<&) denotes a suitable local self-interaction. The transformation from (z,t) to proper time r and rapidity r\ 
is defined as 



t = t cosh 77, t = \Jt 2 — z 2 , dt = cosh 77 d T — - g 

r 

/ z\ . cosh 77 
z = rsmn?7, rj = atanh I - j , o z = — smli 77 o T -\ a^. (2) 



The boundaries in the forward light cone are (note: dzdt — rdrdr]) 



A ft pA ,-atanh^/l-(-r/A) 2 

dt dz=> rdr / dt). (3) 

J -t JO J-atanh\/l-(T/A) 2 

Note that (only) for r = the 77-integral extends from —00 to +00. We therefore consider instead the action 

dr J Cdrj^ J (T-L 2 - U)dr (4) 
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with 

T(r) = ^ (d^fd Vl L 2 (t) = — {d^fdr,, U(t)=t V($)dr,. (5) 

^ J — oo ^ T J — oo J — oo 

Variation with respect to <I> then leads to the equation of motion (EOM) 8S/S& = 0: 

-d T § + <9 TT $ - -^<9„ n $ + 6V/6& = 0. (6) 

T 7~ 

Consider the total energy E(t) at a given proper time r 

E(t) = J {{8S/5d T $) ■ 9 T $ - C)drj = (T(r) + L 2 (r) + U(t)). (7) 

The total derivative of this energy with respect to proper time r is 

dE dE f ( 5E Q _ 8E . 



i(T - £ 2 + 17) + r J (d TT <5> - ^0„„$ + 6V/6$j d T <S> d v . 



(8) 



Using EOM ® leads to 



This is solved by 



with constant C. 



dE I, I, 

— = -(-T-L 2 + U) = —{E-2U). (9) 

CLT T T 



E= — + - I Udr ( 1.0) 

T T 



1. Two remarks: 

(A) : If the potential V(<&) averaged over rapidity is independent of proper time r (as can be realized for random 
distributions of in rapidity space, i.e. during the early phase of an evolution) such that we have U oc r, then a 
simple solution to @ is obtained 

T + L 2 = ^, i.e. £=- + dr (11) 

T T 

with constants Co and C\. The first term ( solving the homogenous part of eq.l©) looks like damping, but it is a 
purely geometrical decrease. (Its time integral is still divergent). 

(B) : If the energy for large proper times r — > oo approaches a finite constant E(t) — > Erx,, then eq.© requires that 

lrnv^oo (T(r) + L 2 (t) - U{r)) = + fluct. (12) 

with eventual fluctuations in U and T + L 2 (which, however have to add up to a smooth total energy E(t), because 
of 4^ — > 0). For proper-time averages over at least one period of the fluctuations we have 

linv^oo < (T(r) + L 2 (t) >= lim^oo < U{t) >= E^/2. (13) 

(Bl): This may be realized such that for large r all time derivatives disappear (T — > 0) and L 2 = U — ► Eoo/2. This 
represents an ensemble of adiabatically shrinking solitons which carry all of the conserved energy E^, (see below). 

(B2): Another way to arrive at an asymptotically constant E appears if the gradient terms L 2 become negligible 
after a long time. Then ea. l|13|l requires that asymptotically < U > and < T > are equal and constant in proper 
time. This implies that there exist non- vanishing average time-derivatives in the system which decrease like I/a/'" in- 
proper time. 
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B. Shrinking solitons at rest in the Bjorken frame 

For a stability condition we omit all time-derivative terms, i.e. put T — and keep in the EOM only 

- -^d m $ + 8V/5<S> = 0. (14) 

This equation describes extended objects whose size in rapidity space shrinks proportional to increasing proper time, 
which here simply acts as a scaling parameter. The derivative of their total energy E = Li + U with respect to proper 
time is given by 

< 15 > 

From multiplying EOM (12) with we have 

1 



9^--(a,^ + y($)j=o, (i6) 

or, if the integrals are finite, 

L 2 = U, and dE/dr = 0. (17) 

This result, that for the shrinking solitons L2 and U are independent of proper time, can of course be obtained by 
simply replacing $(77) by $r(f?) = &(tt]) in the definitions JSJl. 

However, this only holds in an adiabatic limit, where proper time can be considered as a parameter. In fact, for 
the shrinking solitons we have d T i> T = (^/r)9 r; $ r , so 

T=^ Jrf(drMv)) 2 dv. (18) 

So, in the limit of large proper times, T becomes rapidly negligible as compared to the constants Li and U, and E is 
asymptotically conserved. 

C. Including a dissipative term 

Genuine dissipation is taken into acount by adding a dissipative term to the EOM. With dissipation constant 7, 
eq.(6) now reads 

(7 + -)fl!r$ + 3 rr $ - iftw* + SV/5& = 0. (19) 

T T 

Inserting this into eq.JHJl we have 

^ = -~(E - 2U) --it I {d T <5>) 2 d V . (20) 

or t J 

The first term again is purely geometrical, while the rate of energy loss due to dissipation is given by the second term 
alone 

^) = -yr f (d T <S>) 2 d V . (21) 

dT J diss J 

It is equal to —27 times the kinetic energy T, so we have 

d{T + L 2 + U ) = _l_ _ r 

UT T 

For the early stages of an evolution, in analogy to l|ll(l . we expect for the smooth parts of T, Li, and U 

< T + L 2 >= — /(r), <U>=C x t. (23) 
r 



5 



The function /(t) satisfies the differential equation 

% = -2 7 <T> dr. (24) 

The dissipation will exponentially reduce the kinetic gradients as compared to the rapidity gradients, so for the smooth 
parts < T > and < L 2 > we try the ansatz 

< T >= exp(— cryr) < L 2 >, (25) 
where we have allowed for an as yet undetermined constant a. This leads to 

/( T )=(i±f ) (26) 

if we require that /(tq) = 1 for all values of 7, i.e. if for the initial configurations at t = To the integrals T + L 2 are 
independent of 7. This provides an analytical expression for < L 2 > 

C (l + e-°TT)a/«-i 
<^ 2 >________ (27) 

as compared to explicitly solving the EOM (15). We shall see below that its intersection with the linearly rising 
potential < U >= C\t may be used to define a simple estimate for the freeze-out time Tf. 

D. Adiabatic Approximation 

For large values of 7, such that 7x0 ^> 1, all fluctuating terms are suppressed from the very beginning, the evolutions 
get overdamped and we may approximate the EOM (|19|) by omitting all second-order time derivatives. As r > To, it 
is then also sufficient to keep in the dissipative terms only the constant 7 

7 d T $ - -^8^$ + 5V/8& = 0. (28) 



With this equation of motion we now have 



dr t 

Proceeding as in equations 123|) and l|24|l . with U — Cit, we now have 



d{L * + U) -h-L 2 + U)-2 7 T. (29) 



L 2 = ^ exp (-27 y (T/L 2 ) dr 



(30) 



With a suitable ansatz for the ratio T/L 2 during the early part of the evolution we may again obtain a simple 
expression for L 2 and, in connection with U = C\t, for the freeze-out time Tf. 

III. SIMPLE EXAMPLE: $ 4 -MODEL 

The $ 4 -modcl uses the local potential 

n$)^($ 2 -l) 2 , (31) 

so the action S in (z, ^-coordinates is 

S = J Q(^$) 2 - \(d z <S>) 2 - ^($ 2 - l) 2 ) dzdt. (32) 
For stable static solitons we have [2^ 

\ J (d z *fdz = -J{&- l) 2 dz = ^ . (33) 
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\ f (d^fdrj = r 2 \ f ($ 2 - l) 2 ^ = . (34) 



A. Shrinking solitons in the Bjorken frame 

After transforming to the Bjorken frame the stability condition (|14fl looks like the equation for stable static solutions 
in the rest frame with A replaced by At 2 . So we have according to eq. I|33l) 



With definitions ©, we find 

\/2A 

L 2 = U = E/2 = ^— = const T . (35) 
o 

So, in the (adiabatic, i.e. T — > 0) limit of large proper times the total energy of the shrinking solitons in the Bjorken 
frame approaches a constant which coincides with their static energy in the (z, i) -frame. This represents a realization 
of the above remark (Bl) in section II. 

It should, however, be noted that on a lattice these relations are violated if with increasing proper time the shrinking 
size (Ty/X/2)" 1 of the stable solitons approaches the lattice constant (set equal to one). Gradients for single solitons 
then cannot exceed the value of 2 (corresponding to a jump in <E> from -1 to +1 for neighbouring lattice points). So 
we have 

2 

Li — > n— for t — > oo, (36) 

T 

where n is the final number of kinks+antikinks surviving for large proper times. This means that L 2 will not 
approach a constant, but will decrease like 1/t, while the kinks' contribution to U rapidly goes to zero like r -3 . So, 
on a lattice, dE/dr — cannot be reached in an adiabatic limit for large proper times. Still, dE/dr — > for r — > oo 
can be realized by balancing fluctuational contributions to T and U such that E is constant, which constitutes a 
non-adiabatic realization of the remark (B2) discussed in section II. A. 

B. Evolutions on a lattice 

1. Symmetric potential 

Let us first consider evolutions in the symmetric potential 

V {+) {§) = ^ (($ 2 + l) 2 - 1) , (37) 

which has only one single minimum at $ — 0. (We subtract the constant —1 to maintain V(0) = 0). So, a priori, 
we would not expect any 'phase transition' to occur in this potential. Initial configurations at r = tq <C 1 for 
the evolutions consist of values $i at each lattice point i — 0, ...,N, randomly chosen from a normalized Gaussian 
distribution 



/(*) = -A =f exp(~) (38) 



with mean square deviation $ 2 — a 2 . So, the system looks similar for all rapidity intervals that make up the 77 lattice. 
The integral 



l/ (+) ($)dr/ =S> N J y (+) ($)/($)d$ (39) 
for sufficiently small width a <^ 1 is approximated by 

J $ 2 /($)d$ = N^a 2 . (40) 

So, the initial value of t/W at r = tq is given by N^u 2 tq, while J7' + -'(t) will be oscillating around the time-averaged 
< C/ < - + - ) > given by the straight line 

< 17<+) >= N^a 2 r. (41) 
4 
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For an estimate of the integral J (d v &) 2 dri we implement the first derivatives on the lattice as <&i+i — $i and obtain 

For the second equality we have assumed that the two-point correlation function J^^i+i^i factorizes, i.e. that 
$-values at neighbouring lattice points are uncorrelated. So, the initial value for L2 at r = tq is given by 

2 

£2 (r ) = N— . (43) 

TO 

The initial values for the proper time derivatives d T <&i(ro) are not relevant and we can put them to zero. The equations 
of motion immediately after the onset of the evolution build up T such that the smooth sum T + L2 closely follows 
the line Na 2 /r, while T and L2 separately oscillate around the time averaged values 

<L 2 >=<T>=N— . (44) 
2r 

During this early phase of the evolution the potential is unimportant because it is suppressed by two orders of t as 
compared to L 2 . The proper time value 

t) = 2/A (45) 

marks the intersection of < L2 > with the potential < [/(+) > where begins to influence the motion. This marks 
the end of this first (gradient-dominated) phase which consists of freely propagating fluctuations. 




FIG. 1: Evolution on an N = 100 lattice of a single event according to EOM © in the <l> 4 -model with the symmetric 
potential (A = 0.1). The initial distribution at r = to is Gaussian with a — 0.1. The transition time ry is defined by the 
intersection of the lines Na 2 /(2r) and Na 2 \r/4. 

For times r 3> r/, the gradient terms become irrelevant. If, for large r, we neglect in the EOM © the gradient 
terms altogether, then the fluctuations satisfy the simplified EOM 



-d T 4> + d TT cj) + \cf> = 



(46) 
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which is solved by 



-^cos (VAt - 7T / '4) + 0(t- 1 ), (47) 



i.e., the amplitudes of the fluctuations decrease like 1/y/r. This implies that both, and T, oscillate with constant 
amplitude, and frequency 2\/~X, both adding up to an asymptotically constant = T + U^ + \ while Li rapidly 
decreases like 1/r 2 . 

So, for times r 3> 77, we are left with a system dominated by non-propagating fluctuations, or, in other words, with 
a system of N independently moving oscillators, all fluctuating with the same frequency y/X (the 'mass' at $ = 0) 
and amplitudes decreasing like 1/y/r. Asymptotically for large proper times the total energy is conserved, realizing 
case (B2) discussed above. 

Figure n shows the typical features of an evolution in the -potential (for A = 0.1 and a = 0.1). As long as 
To <C Tf, the actual choice of tq and the initial values of (9 T $i(r = To) are unimportant. So, in contrast to our 
original expectations, the evolution shows very distinctly a transition from the early gradient-dominated phase to a 
system of N independent oscillators dominated by the local potential. In the first phase we observe the virial theorem 
< T >=< L2 >, in the second phase we have < T >=< U >. At the transition time t — Tf, which is independent of 
a and To, < L2 > changes continuously from the r _1 to 77 r behaviour. 

The order of magnitude of 77 can be obtained from the following argument: In the Bjorken frame, the field 
fluctuations carry an effective wave number fc/r and mass y/X. If k were not limited from above, then for any r there 
could be propagating waves with (/c/r) 2 > A. However, if the system is characterized by a momentum cut-off, k < A, 
then for r 2 3> A 2 /A, the system decouples into N independent oscillators with the same frequency yf~X. On the lattice, 
k is limited by the finite lattice constant a: k < n/a. So, (with a put to 1), this provides an estimate for the freeze-out 
time Tf. The statistical average over the initial ensemble then leads to l|45[) . From this consideration we conclude 
that the observed freeze-out is the consequence of a momentum cut-off, which in the present case is imposed by the 
lattice. We should, however, keep in mind that here the lattice constant has a well-defined physical meaning: With 
our choice of the initial ensemble, where values of the field variables <&i at neighbouring lattice sites are uncorrelated, 
the lattice constant can be identified with the initial correlation length £. So, in this case, the freeze-out time should 
rather be written as 

r] = jj,. («) 



2. Potential with broken symmetry 

Let us now consider an evolution subject to the potential V (|31|) . which is characterized by two degenerate minima 
at $ = ±1 and satisfies V(±l) = 0. As before, for sufficiently small width a -C 1 of the initial distribution, the 
potential is unimportant during the early phase of an evolution. So, starting from the same initial distribution as 
described above, the first part of an evolution proceeds exactly as in the case of the symmetric potential y(+) 137|l . 
The smooth part of the integral LA - ) — r J V^'dr) = t J (V(&) — V(0)) drj again is approximated by 

< >= {-)nK 2 t, (49) 
while < L2 > is still given by l|44|l. So, the freeze-out time 77 as determined from 

- < cr<-> >=< L 2 > (50) 

again is given by l|45|) . (Due to the opposite sign of < > as compared to the case of the symmetrical potential 

< >, this condition here implies that even for the highest wave numbers k ~ 1/a supported by the lattice, the 
frequencies become imaginary for r > Tf). 

For t > Tf, the gradients which couple the field variables <I>i at neighbouring lattice sites are suppressed and the 
onset of the potential drives each <E>i towards +1 or -1, depending only on the sign of $, at t ~ r/. So, evidently, 77 is 
also that point in time when the number n of kinks+antikinks finally found in the system for r — > 00 is determined. 
Necessarily, because for r > 77 the $i at neighbouring sites no longer interact, the analytical shape of the stable 

(shrinking) kink is irrelevant, or, in other words, the size R = (t\J A / 2^ of the kinks has shrunk below the lattice 
unit. This provides another interpretation of the freeze-out time which we could have chosen alternatively to define 
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FIG. 2: Evolution on an N = 100 lattice of a single event (with a = 0.1) according to EOM JSJ in the <3> 4 -model with the 
broken-symmetry potential V($) Q31[l (for A = 0.1). The freeze-out time Tf = = 4.47 is defined as in figure Q by the 

intersection of the lines Na 2 /(2r) and Na 2 \T/A. The number n of kinks+antikinks remaining in the final configuration is tied 
to the intersection of the lines 2n/r and AT At/4. This particular event leads to n = 24. The arrow at r = Ti points to the 
end-of-roll-down time as given by eq. |JZD- Its intersection with the line NXt/4 marks the finally available energy Eoo (in the 
approximation discussed in subsection III.E). 



Tf. Clearly, for R < 1, kinks and antikinks no longer overlap, so there is no further annihilation and their total 
number n is conserved. 

Figure [3 shows the typical features of an evolution in the symmetry-broken potential V l|31fl (for A = 0.1), starting 
from an initial Gaussian distribution with a = 0.1. For proper times larger than the freeze-out time, r > rt = \j2j\, 
L2 switches from oscillating around Na 2 /(2t) to oscillating around 2n/r (with n — 24 in this particular example), 
while U begins to deviate from NXt/4 and joins with T into combined oscillations around the constant -Eqo/2. These 
fluctuations consist of the N $i-values oscillating independently around the vacua $ = ±1 with frequency ujq = a/2A 
(the 'mass' in the true vacuum) and amplitudes decreasing like l/y/r, so U and T asymptotically oscillate with 
frequency 2ojq with constant amplitudes. 

We can obtain an estimate for the number n by the following consideration: As given in (|36(l , for large proper times 
t > tj we expect L2 to approach the limit < L2 >= 2n/r which on the lattice comprises the kinks' contribution to 
the energy. For proper times up to Tf, t < Tf, the integral J Vdij is dominated by the constant contribution V(0), 
i.e. we have < U >= NXt/4. Both constants, V^(0) before the transition, and 2n after the transition, do not enter 
into the EOM which describes the fluctuations around them and which determines the transition time r/. In normal 
time t, energy conservation would require both constants to be equal. To translate that into the proper-time frame, 
we have to extrapolate < L2 > backwards to the freeze-out time r/ and consider its intersection with < U > at Tf 

< U >=< L 2 > \r=T f i-e. N—Tf = 2n—. (51) 

4 3 Tf 

So we obtain as an estimate for the number n of kinks+antikinks expected to survive, 

N 

n=-. (52) 

For this consideration the fluctuational contributions to U and to L 2 are unimportant because they are suppressed by 
order 0(er 2 ), (apart from the fact that they determine Tf). The result l|52|) is very simple and independent of A (and 
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of cr and To, if chosen sufficiently small). It may be noted that it is just one half of the purely statistical Kibble limit 
tik = N/2. Numerically, it can be measured as the average kink+antikink number n obtained for a large ensemble of 
events. Averaging over 1000 events we find a slight dependence on To and a: for N — 100, a = 0.01 and To = 0.1, 0.01 
we have n = 25.0(4.7), 27.6(5.0); for a = 0.1 and t = 0.1,0.01 we find n = 26.1(4.7), 28.9(5.1), respectively. (Cf. 
figs |6] and |7|). (The numbers in brackets give the corresponding mean deviations). 

Although the momentum cut-off implemented through the lattice (which is decisive for the freeze-out time and 
instrumental for the occurrence of the transition) can have a genuine physical meaning, it is still unsatisfactory that 
the shrinking kinks' contribution to the energy, due to the finite lattice constant, is given by 2ti/t, and not by their 

true energy l|35l) E = n ■ 2^^ which is independent of r. To correct for that lattice artifact we can replace (|51|l by 

N * = 2n ^ (53) 
4 3 y J 

and, with (|45Jl . obtain the result 

n = -N (54) 
8 

which again is independent of A. We shall see in the following that this is the appropriate procedure, if the freeze-out 
occurs before the soliton size has shrunk to the lattice- unit. 



C. Temperature and initial distributions 



We have seen from the preceding results that the momentum cut-off A = ir/a, here imposed on the system by the 
finite lattice constant a, is crucial for the freeze-out time and for the onset of the phase transition. Although, for the 
initial configurations chosen, the lattice constant does acquire physical meaning as the initial correlation length, it may 
seem very unsatisfactory that the essential physical processes rely on a rather arbitrarily chosen momentum cut-off. 
Therefore it is more natural to introduce a cut-off A 7- through a temperature T characterizing the initial distribution. 
For this purpose we construct the initial configurations $i (i = 1, .., N) on the lattice, as Fourier transforms 

iV 



k=l 

of distributions (k=l,..,N) in momentum space, chosen randomly from a Gaussian deviate fa{*&) with fc-dependent 
width at, 

In other words, we choose a Boltzmann distribution for the average occupation numbers < nk >=< >— <j\ for 
particles with mass \[\ ( i.e., the 'mass' of the fluctuations in the symmetric potential ) and temperature T. 

On a lattice, the upper limit for k is N. So, as long as 

2/2 

T < A r , with A T = (56) 

2v A 

the lattice cut-off is unimportant. With a = 1, A = 0.1 we have T < 15.6. Within that limit we obtain for the width 
of the initial distribution of the field configurations <f>j on the 'coordinate' lattice 

(T »=<$>>=Ir<^>=i 1 (57) 

N ^ k 2tt k j 

k 

which even for <7q = 1 still is smaller than 1. So, our previous arguments which lead to | < > \ = Na 2 \/4r still 
hold and we obtain in place of (|49|l 

2 

| < £/(±) > | = NX ^t. (58) 
on 
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For an estimate of < L2 >, as in equation l|42|) we need the correlation function < >, which now no longer 

factorizes, but instead is given by 

- < >= g exp (- VAT/2) . (59) 

i— 1 

So, together with l|57(l we obtain for < L2 > 

< L 2 >= (l - cxp (- VAT/2)) . (60) 
We again use | < [7 < - ± ' > | =< L2 > to get the freeze-out time 

t) = I (l - exp (- VAT/2)) . (61) 
If we use as an estimate for the average kink+antikink multiplicities n the relation (|51(l we find 

n= j (l-exp(-VAT/2)) . (62) 



However, for small temperatures T, the freeze-out time (|61H will be much shorter than the time y/2/X when the kinks' 
size approaches the lattice unit. So it appears more appropriate to replace (|51|) by (|53p. which accounts for the fact 
that at times r < \/2/X the emerging (and shrinking) kinks' energy is constant and given by 2-^^. This leads to 




0.01 0.1 1 10 T 



FIG. 3: Kink+antikink multiplicities n (averaged over 500 events) as function of the initial temperature T (for TV = 100, A = 
0.1). The solid lines are calculated from the expressions 162H and 1631 1 . respectively. 

In figure |21 we compare the numerically measured n (averaged over 500 events, N = 100, A = 0.1) as function of 
initial temperature T with both expressions, (|62|l and Evidently, using the shrinking kink's true constant energy 

2-^p, provides a good approximation for n in the low-temperature region, while the use of the lattice expression 
2/t describes only the limit where the temperature exceeds the cut-off At in (|56|l . These results show that our 
expressions for the freeze-out time and the multiplicities n work reasonably well, independently of the physical origin 
of the momentum cut-off. We shall therefore in the following return to the simple case where the momentum-cut-off 
is provided by the lattice constant. 
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D. Cooling times 



Up to now we have considered the 'sudden quench', i.e. the initial 'hot' configuration was exposed from the outset 
r > tq to a fixed potential. In general, however, the $ 4 -potential would be of a form like 

vm = ^ 2 -v 2 (r)) 2 (64) 

with a time- (or temperature-) dependent function v 2 (t) to drive the phase transition. Typically, the functional 
dependence of v 2 on proper time r would be like 

v 2 {t) = tanh(r/r c - 1), (65) 

parametrized by a cooling time r c which characterizes the time it takes the potential to approach its 'cold' form. 
A possible scenario how the time-scale of this change in the effective chiral potential may be driven by a confining 
transition coupled to it has been discussed in the rapidly expanding Bjorken metric in [21| . In the present context, 
however, we shall simply impose a suitable time dependence like (|65|) . 
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0.01 0.1 1 10 x c 100 

FIG. 4: Kink+antikink multiplicities n (averaged over 500 events) as function of cooling time r c (for N = 100, A = 0.1, and 
for two different choices of the initial width a = 0.0001 and a = 0.01). The frequency of the oscillations around n — N/2 is 
u = \/2A. 



We have seen that in the Bjorken frame for a fixed potential the kinetic freeze-out time r/ = \/2/A in Ij45(l is 
obtained independently of the sign of v 2 , being solely determined by the coupling constant A and the momentum 
cut-off which characterizes the system. Retracing our arguments which have lead to equations (|41() and (|49|l , now for 
the r-dependent potential (j64|) . we obtain for the freeze-out time 

•? - x^y < 66 » 

For an ansatz like l|65|l the solutions t 2 of this equation satisfy r 2 > 2/ A. They approach the lower limit 2/ A for 

t c — > 0. For very long cooling times r c — > oo, we have t//t c — > 1. In any case, if r c is smaller or comparable to y/2/X, 

then also the freeze-out time r/ is comparable to \Jlj\ and we can apply equation Ij51|l which determines the kink 
multiplicities. With the potential l|64|) this is now modified as 

n\v^t s = \{2v 2 ) 2 K (67) 
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So, with ||HSJ|, we obtain the same result n = N/A as before, independently of the value of v 2 {tj). We conclude, that 
as long as the cooling time t c is smaller than the kinetic freeze-out limit \j2j\ the resulting multiplicities are almost 
independent of r c , and it is justified to use the sudden quench approximation r c = 0. 

This is no longer true for r c ^> \/2/X. In that case the potential l|64|l is still symmetric with a single ($ = 0) 
minimum at the time r = \j2j\ when all modes decouple. Therefore, by the time r = r c when the potential 
develops its degenerate true minima with broken symmetry, all field variables <I>i roll down into the nearest minimum, 
independently of each other. So, for large cooling times, the average number n will approach the Kibble limit N/2. 
Figure 0| shows kink+antikink multiplicities n (averaged over 500 events) as function of cooling time t c (for A = 0.1, 
and for two different choices of the initial width a = 0.0001 and a — 0.01). The fact that n oscillates with u> = a/2A 
(the 'cold' mass) shows that the roll down proceeds only after the potential has essentially reached its final 'cold' 
form. 

Physically, th is is not interesting because the field variables $i at different lattice sites do not interact with each 
other for r ^> \/2/\. We therefore in the following consider only the sudden quench r c = 0. 



E. Evolution after freeze-out 



At freeze-out the energy present in the rapidity interval is completely converted into kink-antikink pairs. After 
freeze-out the roll-down sets in, i.e. the field evolution is dominated by the potential. Locally, the field picks up 
kinetic energy T under the influence of the potential — (for |$| -C 1). As long as |$| is sufficiently small the 
equation of motion is 

-$ + $ - A$ = 0, (68) 

T 

solved by 

$(r) = c/ (VAt) = c exp(v / Ar)F(v / Ar). (69) 

with F{x) a slowly monotonically decreasing function of x. If a is the (average) length of <!> at r/, then the constant c 
is c = ct//o(VAtj). For the typical roll-down time t\ we take the time when |$(ti)| = 1 is reached. So, we have 

^ ( n-r /) = -ln. + ln^^l. (70) 

F(VAti) 



As a first approximation we use 



T i ~Tf = -4= m cr- (71) 



(This is a lower limit, because the omitted last term in Q70(l is > 0, and because for |4>| ~ 1 the 4-th order terms in V 
become important. But for an estimate, l|71|l is good enough. The estimate should be best for small a.) During this 
time interval the spatial volume covered by the Bjorken slab grows linearly with proper time. 

The local energy density is dominated by the constant A/4, so by the end of the roll-down the additional energy 
(ti — Tf)NX/A is available to be converted into a- fluctuations around the true vacua $ = ±1. These fluctuations 
carry the mass loq — V2X. So, at the end of the evolution, in addition to the n (anti)kinks formed at freeze-out time, 
we find cr-mesons emitted from the considered rapidity slab with multiplicity 

iVA(-l) 1 Nina 

n a = — —lna—= = =-. (72) 

4 VA V2X iV2 

Altogether, with n = N/A from (|52|l we obtain for the multiplicity ratio of baryon-antibaryon pairs to cr-mesons 

n/2 _ 1 



(73) 



which is independent of the coupling constant A. 

In fig.©, in addition to the arrow at freeze-out time t/ = \j2j\ = 4.47 another arrow shows the end-of-roll-down 
time T\ — Tf + VlOlnlO = 11.75, as given in (|71|) for A = 0.1, a = 0.1. Its intersection with the linearly rising 
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[/(<!> = 0) defines the energy NXti/4 available at r = t\. It can be seen in fig- that the straight horizontal line 
through this intersection coincides very closely with the actual numerically calculated total energy E for times r > t\. 
(The contribution 2jit~ 1 of the kinks is negligible (on the lattice) for r > ri). 

The last equation l|73[) holds for the initial configurations as described in (|38[) . If we use the initial conditions H55|) 
which involve a temperature cut-off in momentum space, we have instead, with l|63[) . (|72|l and (|57l) 

3Vl : (l-exp(-VAT/2)) 1 /2. (74) 



For small values of vAT this reduces to 



41n< 



2/?=- 3 (VXt)1/3i (75) 
n„ 41ner 



With the potential of the true vacuum at $ = ±1 put to zero, after the roll-down there is no further linear rise in 
the total energy of the expanding Bjorken slab. Remarkably, the r _1//2 law l|47|) of the amplitudes of the remaining 
vacuum fluctuations provides an asymptotically constant total energy E(t — > oo) = i?oo 7 similar to the asymptotically 
constant (exact) energy stored in the shrinking kinks. 



F. Including an explicitly dissipative term 

We have seen that in the Bjorken frame the onset of a phase transition, together with the freeze-out of independent 
kinks, antikinks, and 'mesonic' fluctuations around the true vacua, is basically governed by the momentum cut-off 
which characterizes the system. The 'cooling' of the system originates in the decrease of effective wave numbers fc/r 
with increasing proper time r. In addition to this kinetic cooling we normally would expect also genuine physical 
cooling to occur due to radiative processes which dissipate energy from the excited system into the surrounding 
vacuum. 




0.1 1 10 100 x 1000 



FIG. 5: Evolution (N = 100, A = 0.1, To = 0.1, a = 0.1) according to EOM (E3 with additional dissipation (7 = 2). The thick 
solid line shows the expression 1271 (for a = 1/2) for the smooth part < Li >. It intersects | < > | = (jV /4)\o 2 t at 

freeze-out time 77 = 2.05. This event produces n — 16 kinks+antikinks. The constant 16 e with e = 2^^- intersects the total 
energy E near freeze-out time 77. The line 16 • 2/r dominates the energy for r > 200. 
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We incorporate this effect by including a finite dissipation constant 7 into the EOM (|19f) as described in section 
II. C. During the early phase of an evolution, up to the onset of a transition, this does not affect the linear rise of 
< U >. On the other hand, the dissipation damps away the rapidity and proper time gradients such that their smooth 
sum T + L2 drops below the Na 2 /r line. Therefore the freeze-out time 77 is reduced as compared to the case 7 = 0. 
Correspondingly, the average number n of kinks+antikinks produced is also reduced. As an example, figure [5] shows 
the evolution of one event for 7 = 2 (with a = 0.1, to = 0.1). This event leads to n = 16 kinks+antikinks. Applying 
the approximation l|27|) (with a = 1/2) for an average < L2 >, (which in figure [S] is shown by the solid line), we find 
as in l|5U|l from its intersection with the linearly rising | < > \ — (N / '4)A<t 2 t the freeze-out time 77 = 2.05 (as 

compared to 77 = y/2/X — 4.47 for 7 = 0). At this early time the emerging kinks have not yet shrunk to lattice-unit 

size, so we use their true constant energy e = 2-^^ to extract from the total energy E at time 77 their multiplicity 
n = E(rf)/e. In figure |S] this is indicated by the horizontal line at 16 e, which may be seen to intersect E (which 
near 77 is dominated by the linear rise of potential < U >) near proper time 77. The roll-down phase which follows 
for t > 77 is characterized by a peak in the proper time derivatives T, followed by a rapid decrease in U, as the 
configuration settles in the true vacua. Finally the total energy is dominated by the (lattice-artifact kink energies) 
n ■ 2/t behaviour of L2 for n = 16. During and after the roll-down phase all fluctuations have been damped away by 
the dissipative term. 




10 1 1 1 — ^ 1 1 — ^ 1 1 — ^ 1 1 — ^ 1 1 — ^ 1 1 — J 

0.01 0.1 1 10 100 1000 y 10000 

FIG. 6: Numerically measured kink+antikink multiplicity n (averaged over 1000 random events) as functions of damping 
constant 7 for initial time to = 0.1 (and two values for the width of the initial distribution a — 0.1 and a — 0.01) . The solid 
lines show the results (for a = 1/2) of the analytical approximations eqs. 1771 . 1781 and 179H . 180H . respectively. 

This pattern changes for very large values of the dissipation constant 7. As jtq approaches unity, the proper time 
derivatives T get strongly suppressed already from the beginning of the evolution such that T stays small and L2 
dominates the total energy. For 770 3> 1, L2 is determined by its initial value alone, it approaches the line 

L 2 =N— , (76) 

r 

i.e. it takes twice the value l|44|l it has for 7 = 0. Correspondingly, the average number n of kinks+antikinks 
approaches the Kibble limit n — > N/2. In other words, the field configuration is basically frozen from the very 
beginning at t = To, the evolution consists only of the field variables $i at each lattice point drifting to the nearest 
true minimum at <J> — ±1. So, in that limit, the notion of a freeze-out time becomes meaningless. 
Still it is interesting to employ the analytical expression l|27|) for < L2 > to determine 77 from 
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and the number n of kinks+antikinks from i|53|) 



Hr^ (78) 



This provides the limits n — > 3N/8 for 7 — * 0, and n — > 37V\/2/8 for 7x0 1 and a minimum of n near 7T0 ~ 0.5. 
Figure shows the resulting n as function of the dissipation constant 7, for a — 1/2, in comparison with the 
numerically measured n (averaged over 1000 events, for A = 0.1, tq — 0.1 and two values for the width of the initial 
distribution a = 0.1 and a = 0.01). Although the general features of n are nicely reproduced, the limits for 7^0 
and 7T0 3> 1 are not correct. This is no surprise because in both limits the true kink energy e should be replaced by 
the lattice result, i.e. (|78|) should be replaced by ijoTT) 

n = N^rj. (79) 

This then leads to the correct (lattice) limits n — > 7V/4 for 7^0, and n — > N/2 for 7T0 ^> 1. It is interesting to note 
that a quite satisfactory parametrization of the results can be obtained by using (|79(l for all values of 7, but replacing 
< L 2 > on the left-hand side of J77|) by 

1 7 p-°"fTf \ 

-r/. (80) 



77 (1 + e- Q T T ») 2 



This comparison is also shown in figEl again for a = 1/2. Evidently, it provides remarkably good agreement for 
7- values up to 710 < 5. 



G. Adiabatic approximation 



For large values of the dissipation constant 7 it should be sufficient to employ the adiabatic approximation as 
described in section II. D. Evolving the initial configurations with the purely dissipative EOM (|28|l . leads to the 
kink+antikink multiplicities n (averaged over 1000 events, for TV = 100, and A = 0.1) shown in figuredfor two initial 
times To = 0.1,0.01 (and initial Gaussian width a = 0.01). For comparison, the non-adiabatic results for the same 
two cases are included. 




FIG. 7: Numerically measured kink+antikink multiplicities n (averaged over 1000 random events) as functions of damping 
constant 7 for two different initial times to — 0.1, 0.01 ( and initial width a = 0.01 ) as obtained from the adiabatic EOM 1281 . 
For comparison, the non-adiabatic results are also included. The solid lines show the result of the analytical approximation 
1791 . ISOI (with a = 1/2) for the same two cases. 
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It is apparent that the adiabatic approximation works very well for jtq > 1. It can also be observed that for 
jt < 1CP 2 the adiabatic evolutions produce no kinks at all, i.e. the configurations become so smooth at a very early 
time (when the 'size' of the kinks is still very large), such that the complete configuration rolls down into the same 
true minimum. The solid lines show the result of the analytical approximation (|79|l . (|80l) for both cases To = 0.1,0.01. 

IV. CONCLUSION 

The dynamics of interacting fields in rapidly expanding systems is characterized by effective wave numbers decreasing 
with time. This is conveniently modeled by transforming from space-time (x,t) to the rapidity - proper-time (77, r) 
(Bjorken) frame. In the equations of motion for the evolving fields, the rapidity gradients then explicitly carry 
the inverse proper-time factors, and an additional first-order proper-time derivative term appears which resembles 
dissipative dynamics, although no genuine physical dissipation is present. In the (?J,t) frame, the system still is 
characterized by a lagrangian which, however, explicitly depends on proper time. So, the energy is not conserved, but 
may decrease or rise during the evolution of a field configuration. This is especially interesting for those cases where 
the system is able to undergo a 'phase transition', where initial configurations through a cooling process settle down 
into the true vacua. Usually, cooling or heating is provided by energy exchange with external degrees of freedom (heat 
bath), but in the (77, t) frame the transition is solely driven by proper time. 

Systems characterized by multiple different degenerate vacua are of specific interest because during the course of the 
phase transition it must be decided how the different vacua get populated. Field configurations connecting different 
vacua will finally show up as defects, kinks, or solitons, and their multiplicities can be physically significant and 
observable. 

The mechanism of how the transition is effectuated by proper time in the Bjorken frame is quite transparent. The 
term which comprises the (rapidity-)gradients contribution Li to the total energy carries a factor 1/r, while the 
potential-energy part U carries a factor r. Therefore, the evolutions at very early times will be dominated by Li 
independently of the potential, while, if the gradients in the field configurations are limited by a momentum-cut-off A, 
then for very late times the gradients will become unimportant, propagation in rapidity space will end, and the field 
variables evolve locally according to the local potential. Evidently, there will exist a freeze-out time 77 which separates 
both regimes. It is conveniently defined as the time when Li equals U (where suitable care has to be taken to subtract 
eventual constants in the definition of the interaction). Clearly, the momentum cut-off ( or the finiteness of Li ) is 
a necessary condition for the transition to occur, and the actual value of the freeze-out time will depend decisively 
on A. Therefore it is important that the momentum cut-off reflects a genuine physical characteristic of the initial 
field ensemble, typically an initial temperature or inverse correlation length. This then, however, implies that the 
freeze-out time can be calculated very simply from the initial ensemble. 

Stable solitons in the Bjorken frame shrink with increasing proper time, or, in other words, their stable size relative 
to the expanding system decreases. But remarkably, their total energy is constant and equal to their static total energy 
in the usual (z,t) frame. This allows to estimate their multiplicities from the total energy available at freeze-out time 
(including eventual constants in the definition of the interaction). This means that we may obtain estimates for the 
average soliton(+antisoliton) multiplicities directly from the form of the initial ensemble and the static lagrangian 
which determines their mass. 

A convenient way to introduce a momentum cut-off is to implement the system on a lattice. Then, however, we 
have to make sure that the lattice cut-off does not interfere with the increasing gradients which characterize the 
shrinking solitons at t/. Alternatively, we have to use the lattice expression for the soliton's energy if it has shrunk 
to lattice-unit size. 

We have investigated these features in the (l+l)-dimensional $ 4 -model. We find that the simple analytical ex- 
pressions for Tf and n agree very nicely with numerically measured freeze-out times and kink+antikink multiplicities. 
We have also used Boltzmann ensembles with temperatures well below and up to the lattice cut-off. Employing the 
appropriate correlation function we again find good agreement with the numerical averages. For low temperatures we 
obtain a square-root law n cx vf. For temperatures approaching the lattice cut-off the influence of the appropriate 
soliton energy becomes evident. We show that within the physically interesting region it is sufficient to consider 
the sudden quench where the initial ensemble is exposed from the beginning to the fixed cold potential. After the 
freeze-out the additional energy picked up by the expanding rapidity slice is converted into mesonic fluctuations. A 
simple estimate of the roll-down time provides an expression for their average multiplicities. 

In addition to the 'free' evolutions in the Bjorken frame we have also investigated the inclusion of additional genuine 
dissipation with variable damping constant 7. It is very satisfactory to find that all the previous considerations can 
be applied to this case, that the decrease in the freeze-out times caused by the damping can be represented in simple 
analytical expressions and that the multiplicities n again reflect the soliton energies at freeze-out time. Finally we 
find that the validity of the adiabatic approximation, (where all second order time-derivatives are dropped), is limited 



18 



by 7T0 > 1 where To denotes the proper time at which the initial ensemble is prepared and the evolution starts. 

It will be most interesting to look at the corresponding features in the (3+l)-dimcnsional 0(4)-model considered 
previously as a model for baryon-antibaryon production in relativistic heavy-ion collisions and compare the results 
with recent experimentally observed nucleon-antinucleon multiplicities and with statistical models. The scheme may 
also be useful for applications in cosmology. 
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